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Abstract. The dynamics of TV > 3 interacting particles is investigated in the 
non-relativistic context of the Barbour-Bertotti theories. The reduction process 
on this constrained system yields a Lagrangian in the form of a Riemannian line 
element. The involved metric, degenerate in the flat configuration space, is the first 
fundamental form of the space of orbits of translations and rotations (the Leibniz 
group) . The Riemann tensor and the scalar curvature are computed by a generalized 
Gauss formula in terms of the vorticity tensors of generators of the rotations. The 
curvature scalar is further given in terms of the principal moments of inertia of the 
system. Line configurations are singular for TV ^ 3. A comparison with similar 
methods in molecular dynamics is traced. 



1. Introduction 

Since the very gestation of the Newtonian mechanics Q , severe criticism of its concepts 
of absolute space and time has proliferated. The criticism concerning absolute space, 
as the tool determining the privileged inertial frames is best expressed by Mach's 
principle ||: inertial frames are an imprint of the matter content of the universe. 
According to Leibniz ||, time is nothing more than the sequence of events. Some 
of these ideas were incorporated in the general theory of relativity. There, time is 
not absolute and Mach's principle is satisfied at least for closed universes (| || . The 
generic idea of timelessness of both a non-relativistic and a general relativistic complete 
universe was analyzed recently by Barbour [||. 

In this generic framework there were also several attempts to incorporate Mach's 
principle and Leibniz's concept of time in dynamical theories at the non-relativistic 
level. Newtonian absolute space is eliminated from such models and in some of them 
absolute time too, by implementing time reparametrization invariance. However, in 
contrast with general relativity, these models retain absulute simultaneity. The study 
of such models dates back to the early attempts of Reifiner || , who has considered 
the dynamics of the system of N interacting particles which depend on N(N — l)/2 
mutual distances. 

In the model developed by Lynden-Bell both the kinetic and the potential 
term in the Lagrangian are invariant under time- dependent translations and rotations. 
This theory, as described by Lynden-Bell and Katz |I| is equivalent with Newtonian 
mechanics in a universe of zero net angular momentum about its barycenter, provided 
the latter is in uniform motion. There exists a strong evidence for the assumption 
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of vanishing angular momentum of the universe, given by the analysis of the cosmic 
microwave background radiation [ flO| , |Tl| |l2| , in support to the validity of the theory. 
Lynden-Bell's model describes an iV— body system of interacting particles in terms of 
purely relative quantities. No overall (time-dependent) rotation of the frame about 
the barycenter or translational motion of the barycenter has any effect on the action 
functional or the equations of motion derived from it. This happens because no other 
distances and velocities appear in the action than the distances between pairs of 
particles and the differences of angular velocities of the lines joining these pairs of 
particles, respectively. General relativistic extensions of these ideas were discussed by 
Lynden-Bell, Katz and Bicak ||. In Lynden-Bell's model, however, time is absolute. 

The fully relational dynamical models introduced by Barbour and Bertotti Q 
are more general than Lynden-Bell's theory in bearing invariance with respect to the 
whole Leibniz group (including a monotonic redefinition of the time): 

x a ' = K a ' b (t)x b + A a '{t) 

t'=t'(t). (1) 

Here 1Z is an orthogonal matrix and A a are arbitrary translation parameters. As a 
consequence of employing an arbitrary time parameter, the Lagrangian for a system 
of N particles becomes a homogeneous function of degree one in the "velocities" 
(dq.i/dt' , dt/dt'), implying via Euler's theorem a vanishing Hamiltonian for the system 
[PI . As an immediate consequence of this and the invariance (1) seven constants of 
the motion identically vanish. These are the energy, the linear and angular momenta. 
Barbour and Bertotti start from "the relative motion of the universe treated as a single 
entity and then recover the motion of subsystems within the background provided by 
the universe at large" .The particular product Lagrangian proposed in [ p"3| although 
being able to (qualitatively) reproduce some general relativistic effects, like the 
perihelion shift and the Lense-Thirring effect, leads also to unpleasant consequences: 
anisotropic masses (in contradiction with experiment), time-dependent gravitational 
constant and a violation of the Birkhoff theorem. In their second paper, Barbour and 
Bertotti |l4| have proposed a more generic framework for deriving intrinsic dynamics, 
relying on the concept of intrinsic differential. Both Galilei- and Lorentz-invariancc 
could be derived for particular subsystems of the whole Leibniz-invariant universe. 

As the theories of Barbour and Bertotti share with the theory of general relativity 
the property of possessing both the Hamiltonian and momenta constraints, the study 
of these models, particularly the issues concerning quantization, can give an insight of 
how such theories should be handled in a quantization process. Particularly interesting 
are the solutions given to conceptual and interpretational difficulties by Barbour and 



Smolin 16 , where an explicit solution is presented for the case of three particles in 
one dimension. In a second paper jl7|], Barbour and Smolin introduce the concept 
of variety in order to provide a (non-local) description of dynamical systems without 
any reference to elements lying outside the domain of valability of these theories. 
The Barbour-Bertotti theories are among the toy-models on which the applicability 
of their new concept and principle are checked. A group-theoretical quantization of 
the Barbour-Bertotti theories was done by Rovelli jl8|], eliminating any concept of 
time. In contrast, an approach relying on a quantization in terms of an intrinsic time 
was performed by Gambini and Mora |t9|. Arguing that time is a manifestation of 
dynamics, Barbour arrives to the notion of timelessness of quantum gravity Q and 
presents an approach to recover time from a static wavefunction of the universe ]2C| . 
We will not enter in the details of the quantization process here. Our aim in this 
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and the subsequent companion paper |2jj is to analyze several aspects related to the 
geometric structure of these theories. 

In what follows, by a Barbour-Bertotti type theory we mean a nonrelativistic 
dynamical theory obeying the Hamiltonian, momenta and angular momenta 
constraints. In Sec. 2 we start from the Newtonian description of the 
interacting N— particle system, imposing the above constraints. Then, by solving 
the constraints and eliminating the Lagrange-multiplicrs wc obtain the reduced 
Lagrangian, homogeneous in the velocities. 

This action determines a Riemannian metric, in a form of a line element of some 
Riemannian space. Extremizing the action is equivalent with finding the geodesic 
motions. These issues are discussed in Sec. 3. Although the particular form of the 
Lagrangian differs from that given inJl3J, it has a product structure, in accordance 
with the general framework settled in |l3| , 14[ . What we recover is the Jacobi principle 
for the constrained N~ particle system, as expected from the argumentation in g. 
General relativity, in the form given by Baierlein, Sharp and Wheeler |p2f resembles 
also a Jacobi principle rather than a parametrized particle dynamics [|| . 

In Sec. 4. we study the geometry of the reduced space for the generic case of 
N > 3 particles. (For N = 1 there is no relative distance at all. For N = 2, as the 
Lagrangian is homogeneous of degree one in the single relative velocity, the Euler- 
Lagrange equation is trivially satisfied and there is again no relative dynamics |Tq] .) 
For N > 3 we find that the metric is degenerate and the degeneracy directions are 
related to the generators of translations and of rotations in the physical space. This 
is again in accordance with |l3| , p"4| , where dynamics is defined on the space of orbits. 
We compute the curvature scalar of the 3N — 6 dimensional reduced space (the space 
of orbits) by means of a generalized Gauss equation, deduced here. The curvature 
scalar is expressed in terms of the principal moments of inertia of the system. This 
allows us to conclude that the line configurations are curvature singularities for N =/= 3. 
We will investigate the N = 3 case in detail in the companion paper pj[ . 

There are certain connections of our work with the analysis of the N— body 
problem by geometric methods in molecular dynamics. A recent summary of this 
topic was given by Littlejohn and Rcinsch (25). We trace a comparison with their 
results in Sec. 5. 

There are several type of indices appearing in our formulae. The coordinates of the 
configuration space are labeled by a particle index, denoted by a capital letter running 
from 1 to N, and a Cartesian coordinate index of the physical space (abstract index), 
denoted by a Latin letter, which takes the values 1,2,3. We adopt the convention 
to write the summation over particles explicitly and apply the summation convention 
for generic coordinate indices. The (Latin) indices of the Kronecker delta 8 a b and of 
the totally antisymmetric Levi-Civita symbol e a b c are raised when necessary in order 
the summation convention to apply. Kinematical and dynamical quantities written in 
a system of coordinates originating in the center of mass carry a distinguishing zero 
subscript, like x Ac . There is one exception over this rule: the principal moments of 

o 

inertia 1^ are written without the zero subscript. Greek letters label coordinates in 
the system of principal axes of inertia, with the origin in the center of mass. The 
summation over these indices will be also explicitly indicated. 
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S[x,p] 



2. The reduced Lagrangian of the Barbour-Bertotti model 

We start from the canonical form of the action characterizing a system of N interacting 
particles: 

where x Aa , dx Aa /dt' and pAa are the coordinates, velocity and momenta components, 
respectively of the particle with mass tua, and the Hamiltonian is given by: 

N 1 

H = Y P AaPAbS ab + V . (3) 

A=l A 

The potential term V is a superposition of the two-particle potentials Vab, which 
depend only on the relative separation between particles A and B: 

V=yV AB (\v A -r B \) , (4) 

A<B 

r A = {x Aa } being the position vector of particle A. 

The Hamiltonian (|3j) is time-independent, implying a conserved energy. The total 
linear and angular momenta 

N N 

Pa = yPAa , L a = y tabX Ah PAc (5) 

A=l A=l 

are also conserved: 

{P a ,H} = 0, {L a ,H} = 0. (6) 
The above constants of the motion form a closed algebra: 

{P a , P b } = , {L a , L b } = e ab c L c , {P ai L b } = e ab c P c (7) 

The last group of the above Poisson brackets bears the message that the total momenta 
transform as vectors under rotations. 

The Barbour-Bertotti theory is obtained from the above Newtonian theory by 
imposing the constraints! : 

P a = 0, L a = 0, H = 0. (8) 

Among them, the condition P a = is fulfilled by any system if the origin of the frame 
is chosen in the center of mass. The other two constraints do not hold in generic 
situations, for example in the case of an interacting system of harmonic oscillators, 
where the energy is positive definite. For attractive forces H = is the marginally 
bounded situation. 

We add to the canonical action (|^) the linear and angular momenta constraints 
(^) with Lagrange multipliers N a and M a . By replacing the time t' with a labeling 
parameter t pa] , the Hamiltonian also gets a multiplier N: 



S[x, P ;M,M a ,M a ] = J dt (^2p Aa x Aa -AfH-JV a P a -M a L,)j 



(9) 



f Without imposing H = and when Vab ar e the Newtonian gravitational potentials, the Lyndcn- 
Bell theory |9l is recovered. 



(The notation x Aa = dx Aa /dt was introduced.) Variations of this action as a functional 
of the multipliers generate the constraints. Variations of (^|) regarded as functional 
of the coordinates and momenta give the equations of motion, equivalent with the 
Newtonian equations for a system satisfying the constraints. 

Next we proceed with the reduction of the action by solving the constraints and 
expressing the Lagrange multipliers as functions of configuration data. First by varying 
PAa in @ we find the relation between the velocities and linear momenta: 

± Aa = M PAb 6 ab + j^fa _ ^ Ab (10) 

m A 

Here we have introduced the dualized form of A4 a , defined by the relations: 

M ab = e abc M c , M a = \e abc M ab . (11) 



The unusual relation (|1C|) between the velocities and momenta is consequence of the 
constrained nature of the system (j^). From Eq. (|o|) we express the momenta as: 

PAa=N- l m A (i Ab 5 ab -N a +M ab x Ab ) . (12) 

The total linear momentum of the system is found by summing up the above momenta 
over all particles: 

P a =Af- 1 M(x b 5 ab -Af a +M ab x b ) . (13) 

Here M — m A is the total mass of the system, x a = Xm=i TUAX Aa /M and x a 

are the coordinates and the velocities of the center of mass, respectively. 

In a similar manner by summing up over all particles the angular momenta 
l Aa = £ ab x PAc, the total angular momentum emerges: 

L a =N- 1 (l a -I ab M b -MN ab x b ) . (14) 

In the above expression we have introduced the velocity-based total angular 
momentum: 

N 

L = e abc m A x Ab i Ac , (15) 

A=l 



the tensor of inertia: 



I 9 h = e l cg e ldh m A x Ac x Ad (16) 

A=\ 

and the dualized Lagrange-multiplier Af ab related to Af a as: 

N ab = e abc Af c , M a = \e abc N ab . (17) 

By virtue of the constraints (||) the equations (|l^) and ( |l4] ) form an inhomogeneous 
linear algebraic system in the multipliers AT a and Ai a . We transform this system to 
an equivalent form by expressing the multipliers J\f a from ( |l3| ) then inserting them 
into ( [l4] ) obtaining: 

K = M ab x b + x b 5 ab , I ab M b = l a . (18) 



() 



Here 

1 N 

x Aa = x Aa _ ^ = /a ± y Aa (19) 

M ^ y ' 

A=\ 

N 

l a = e abc y / m A x Ab i Ac (20) 

c — ' o o 

A=l 

N 

1 g h=e l cq eidhY^m A x Ac x Ad (21) 

o a ^ — ' o o 

are the coordinates, the velocity-based angular momentum and the inertia tensor in 
the center of mass frame, respectively. 

There is a solution of the second equation ( [l8| ) only when the tensor of inertia in 
the center of mass frame is nondegenerate f . A careful study shows that a degenerate 
inertia tensor would only occur in the case of all N particles being along a line (giving 
a rank 2 inertia tensor) or in the unphysical situation of all particles at the same 
point (rank inertia tensor). The inverse inertia tensor (I ~ 1 ) ab is well defined for all 

o 

other configurations of the particles. Thus, modulo the above exceptional cases, the 
multipliers A4 a are: 

M a = {l- 1 ) ab l b . (22) 



If we insert the expressions of the multipliers J\f a and M a from (|l8|) and (|22j) in the 
expression of the momenta ( |l2| ) we obtain: 

PAa = Af- 1 m A {± Ab S ab + e abc (l - 1 ) cd l d x Ab ). (23) 

o 0° 

We have solved the linear and angular momenta constraints, thus the terms 
containing the multipliers N a and M. a can be dropped from the action By 
inserting the expression of the momenta (p3|), the remaining Liouville form and lapse- 
Hamiltonian term become: 

A / 1 



A=l 

1 



. PAai Aa = 2W T - - 1 a bM a M (24) 

* — \ o Z o 



NH=N- l [T -^-IabM a M b )+NV . (25) 

\o 2 o I 

Here T is the kinetic energy in the center of mass frame: 

o 

1 N 

T = - Y / m A x Aa x M 5 ab . (26) 

A=l 

The action is now a functional on the configuration space and function of the lapse 
alone: 

S[x;Af}= JdtL, L{x,x;Af) =AT 1 (V -^I ab M a M b ^j -AfV .(27) 

f One's restricted ability in expressing the Lagrange multipliers in terms of the configuration data 
arises also in the thin sandwich conjecture of gener al r elativity. The local solvability of that 
problem for a large class of initial data was proved in |23| . Recentl y r esults on the existence of 
solutions of the generalized thin-sandwich problem were also obtained |24| . 
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The variation with respect to Af gives an equation which we use to express the lapse: 



T - \l ab M a M b IT — \ {l~ 1 ) ah I alb 
N =\j° °— = ]J2 °_ ° ° , (28) 

By eliminating all multipliers, the task of reducing the action is completed: 



S[x] = J dtL RED , L RED (x, x) = -V (t-\{I - 1 ) ab l a l o t) • (29) 

We have obtained the Jacobi Principle version of Lyndcn-BelPs purely relative 
Lagrangian with the energy set to zero. 

3. The induced geometry 



The action (|29|) obtained by completing the reduction process (for all noncollinear 
configurations) is homogeneous of degree one in the velocities. Therefore the integrand 
can be regarded as a line element in some Riemannian space. In order to make this 



manifest, we write the reduced Lagrangian (29) in two alternative forms 



L RED {x,x) = J-2V(x) GAaBb(x)± Aa X Bb 

V L o o 

= y/-2V{x) [G A aBb(x)x Aa x Bb ] , (30) 
where G AaBb & n d GAaBb Riemannian metrics denned as: 

o 

GAaBb = m A S A BS a b - m A m B (I - 1 ) ef e eca e fdb x Ac x Bd (31) 

O O O O 

m A m B x , , 

LrAaBb = GAaBb 77 O ao . (6Z) 

Without changing the notations, from now on an index pair Aa is viewed as a single 
index running from 1 to 3iV. Though the reduction process at the Lagrangian level 
was accomplished, the configuration space is still 37V-dimensional. It will be one of 
the tasks of the next section to show that the 3iV-dimensional configuration space is 
superfluous and the (3iV — 6)-dimensional space of orbits of the Leibniz group is better 
suited as a configuration space. 

The metric has the simpler form (31) in the center of mass frame, however this 
expression is covariant only under rotations in the physical space. In contrast, the 
expression (32), containing an additional term, is covariant under both translations 
and rotations in the physical space. It can be brought into a completely generic 
form by inserting the expressions (19) of the coordinates x Ac in (32) and in (j" 1 )^ 

o o 

appearing there. 

The geometry underlying the physical motions is characterized by the conformally 
scaled (Jacobi) metric — 2V(x)GAaBb- Extremizing L RE d is equivalent with looking 
for the geodesic motions associated with this metric. Dynamics arise from geometry. 
To our knowledge it was Synge |26| who first applied this viewpoint in a discussion 
of conservative systems in the framework of a "geometro-dynamical theory of the 
manifold of configurations". A recent discussion of the Jacobi metric for simple 
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dynamical systems [£7j concentrates on the singular regions at E — V (at V = 
when specialized to our case). 

As V(x) varies from case to case, in the remaining part of the paper we investigate 
the Riemannian metric (32). Formally, the metric GAaBb is characterizing the free 
motions pertinent to the V = const case. 

At the end of this section we stress that the elimination of the multiplier Af from 
the action ( p7| ) implies a choice of time. This is because the elimination was carried out 
by employing 6S[x; Af]/8N which is the Hamiltonian constraint H = T + V = E = 0. 
However in a relational theory no time exists a priori, in consequence the energy 
conservation equation is not giving the speed of the system on its trajectory in the 
configuration space, as usually. Rather it defines a unic time t for all subsystems |J. 



4. The reduction to the space of orbits 

For infinitesimal translations and rotations dx Aa = d£ a + t a hc x Ah drf we have 
GAaBbdx Aa = . Therefore the metric is degenerated, the directions of degeneracies 
being the vector flows in the configuration space induced by translations and rotations 
in the physical space. 

With respect to the flat metric 

9AaBb = m A 6 AB 6 ab , g AaBb = —S AB S ab , (33) 

m A 

the generators of translations ^ = S b and of rotations rj Bb = e b C jX Bc do not form an 
orthogonal set: 

A,B=1 

E SAaBbtfavffi = tic,* 1 (34) 

A,B=1 
Af 

E 9AaBbV^)Vu) = h 



A,B = 1 



excepting the case when the coordinate axes originate in the center of mass (x c = 
and Iij — Iij) and they are chosen along the principal axes of the inertia tensor 

o 

(where / M „ = 1^5^, 1^ denoting the principal moments of inertia). In such a system 

o 

an orthonormal set of six generators with respect to the flat metric ( |33"| ) is given by: 
7 Bf3 - — !— rV 3 w BP - —ft 3 r B ^ (^) 



The generators (|35|) correspond to translations along, and rotations about the principal 
axes of inertia, originating in the center of mass. The dual set of these generators is: 

«W = ^, v { ll = ^ a ^. (36) 



The expressions of the generators in generic frames are obtained by the tensor 
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transformation law: 



7 Bb _ 
z (m) - 



E- 



dx 



Bb 



Pi 

o 

dx c ~< 



U Aa ~ U C 7 Q x Aa 

C 7 o 



' J Aa 



Bb 



Cj o 



dx c ~< 



C 7 o 



(37) 



We conclude that GAaBb is the metric in the (3iV — 6) dimensional space of orbits 
of the group of translations and rotations, rigged by the generators (37) f : 



G A aBb = 9AaBb ~ ^ U Aa U Bb ~ 2^ V Aa 



Bb 



(38) 



We define 



*B6 - d B d b - 2^ Z (p) U Bb - 2^ W {v) V Bb 



(39) 



It is an easy exercise to verify that Pgff is a projection operator 



ryAa r>Cc T>Aa 

r Cc r Bb — r Bb 



(40) 



to the subspace rigged by the generators: 



rjAa Bb 
Bb z (p) 



= P, 



Aa n M 



Bb a Aa 



= 



P,Aa Bb 
Bb W („) 



_ pAa \y) _ - 
— -TBb v Aa — U 



(41) 



The construction we have found is quite generic. Whenever one has a Riemannnian 
manifold and a gauge group of isometries acting on it, the metric on the quotient space 
is naturally given by the restriction of the initial metric to the orthogonal complement 
to the gauge group. We stress that none of w?Jl and v^l an d neither GAaBb are well- 
defined when any of the principal moments of inertia vanish. Thus we have to limit 
the validity of the above construction to the noncollinear configurations. The reduced 
space is the quotient space R 3N /E(3) where E(3) = R 3 xi 50(3). By extending 
this manifold to include the collinear configurations either, we find a manifold with 
dimensionality 37V — 6 with boundary for TV = 3 and without boundary for N > 4 
p5[ . The case N = 3 will be discussed in detail in the companion paper plj . For 
N = 4 the manifold in discussion is homeomorphic with R 6 |2^] . 

We denote by V and V the connections compatible with the metric G and the 
flat metric g. The respective covariant derivatives of an arbitrary vector field q Aa in 
the coordinate basis { g ® B b } will be denoted by Vccq Aa and Vccq Aa - For each of the 
generators (53) we define one of the tensors: 



(/■0 

XAaBb 



J AaBb 



N 



E,F=1 



Ee r>Ff 



N 

E 

B,F=1 



Ff 



^Aa^Bb V EeV F f 



(42) 



Pairs of indices Aa of these tensors are raised and lowered with the flat metric (E33|) . 



f For a generic discussion of a rigged (n — m) dimensional manifold in an n— dimensional manifold 
endowed with a linear connection see Schouten 
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From straightforward algebra on Eq. (21) written in the principal axis frame the 
following useful expression emerges: 

N W v 

£ m A x A » x A " = — ^(1 - 25„)I P . (43) 

A=l p 



By use of the above relation, together with Eqs. (|35|), (36) and (19), we compute the 
expressions for the x and tensors in the principal axis frame in terms of the 
principal moments of inertia: 

X-AaB/3 — u 

X») _ m B ^ ( Sab - mA 



a/3 



m A m B s^ x b 7x A6 



M 



7<5p 



v(E^)+^(E 



cr A r N 

E <?\ (E e -^ m ) (E ^) (E m c % Gi1 % Gx 



(44) 



\p X ° T G=l 

The tensor transformation law gives the expression of lj^ v ' in a generic frame: 

dx Ca dx D $ 

(u) X " (u) o o (AK\ 

^AaBb — W CqD/3 q x Aa Q x Bb ' 

CDa/3 o o 

Similarly, X^AaBb = ^ h°lds in any frame. 

The tensor X^AaBb being symmetric (it vanishes) is the extrinsic curvature of 
the hypersurfaces with normal vectors u Aa . By contrast, the tensor to^ Bb is not 
hypersurface orthogonal, being antisymmetric. It is the vorticity tensor for the vectors 
w Aa . This antisymmetry originates in the transformation properties of w B ^ under 
reflections: it behaves as an axial (pseudo) vector rather then a polar (true) vector. 
The antisymmetry of the tensor uJ^aBb complicates the relation between the Riemann 
curvature tensors FL A °BbCcDd (here vanishing) and R Ao BbCcDd-> associated to the metrics 
9AaBb and G Aa Bb, respectively. 

We derive the required generalized Gauss equation as follows. Let q Bb , r and s Dd 
be three arbitrary vectors and p Aa an arbitrary one-form in the tangent and cotangent 
spaces of the space of orbits, respectively. The two Riemann tensors are defined as: 

N 

^Aa 



E R A BbCcDd = 2V[ Cc V_ D(i ]<3 / 



B=l 
N 

E RA BbCcD d I*" = W [Cc V Dd] q Aa . (46) 
B=l 

Evaluating the expression R A °BbCcDd PAa,q Bb r Cc s Dd , we get the desired relation 
between the Riemann curvature tensors of the metric g on the 37V dimensional space 
and of the metric G on the rigged (37V — 6) dimensional space: 

JV 

EpEe pAapFf p Gg p Hh Bb _ M (p Aa\ 

n FfGgHh^Ee^Bb ^Cc^Dd 1 - z £^ ^[CcDd] 1 ) 

B,E,F,G,H=1 v 
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N 

E 



nAa 9 V^ / (v) Ao (v) M (y) Aa 

n BbCcDd L V°\Cc UJ Dd]Bb ^[CcDd^Bb 



q " ,(47) 



where l w , u) q denotes the Lie-derivative C w , u) of the vector q Aa , projected to the 
reduced space: 

£ WM q Aa = P£Z£ WM q Bb . (48) 

The major inconvenience induced by the antisymmetric oj^l Bb tensors in (47) is the 
presence of the derivative term of q Bb . For generic riggings such a term should be 
absorbed in the definition of the Riemann tensor of the reduced space as described 
by Schouten p8f . However in our case the situation is simpler. As translations 
and rotations are symmetries of the system, the quantities defined on the space of 
orbits (in particular the components of the tangent vector q Aa ) should not depend 
on the parameters of translations and rotations. In an adapted coordinate system 
the Lie-derivatives are partial derivatives with respect to these parameters, thus the 
Lie-derivative terms in (47) vanish. Because the vanishing of a tensorial quantity is a 
invariant statement, this result holds in arbitrary coordinate system as well. Therefore 
and by employing the flatness of the metric <?, the expression of the Riemannian 
curvature tensor of the rigged space emerges purely in terms of the vorticity tensors 

pAa _ ( ,0) Aa (u) (v) M Aa\ ,. Q x 

K BbCcDd - Z 2^ V W [Cc ^Dd\Bb~ ^CcDd^Bb ) • \ m ) 
v 

This is the desired generalized Gauss equation. 
The curvature scalar is then readily found: 

JV 

H - G H BbAaDd - ~ 6 Z < Bb 10 Aa • \ bU ) 

A,B,D=1 v 

Here G BbDd denotes the inverse metric of the rigged space. A lengthy, but 
straightforward computation, performed in the system of principal axes of inertia, 
employing the explicit expressions of the vorticity tensors (44) and the auxiliary 
relation ( [43|) has given the scalar curvature in terms of the principal moments of 
inertia and the number of particles: 

s = 6 < w - 2 >i:i-^(i>r (51) 

R being scalar, the above expression is independent of the particular frame employed 
for its evaluation. 

The concise form ( |5l| ) allows one to find out whether the manifold of orbits can be 
extended by the inclusion of the alligned configurations. The metric GAaBb extended 
there by a limiting process, is singular. The question arises, whether these line 
configurations represent coordinate or real singularities. Suppose that the system 
is passing through a sequence of configurations C n toward the line configuration Cu ne 
with all particles aligned on the x-axis. There 1% = while I2 = J3 = I are finite. The 
curvature scalar (|l|) behaves as 

hm fl= 6(2*-5) + 6(*-3) 
c„^c Hne I h 

Excepting the case N = 3 the curvature scalar diverges together with the vanishing 
of I\. Thus the line configurations are curvature singularities for N ^ 3. 
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5. Connections with molecular dynamics 

It is instructive to compare our description with a geometric treatment of the N— body 
problem in the language of fiber bundles, employed in both classical and quantum 
molecular dynamics. A recent review of this topic was given by Littlejohn and Reinsch 
p5[ . In that approach the translational degrees of freedom are subtracted from the 
very beginning by the introduction of the 3iV — 3 Jacobi coordinates. When the 
configurations are not collinear, the remaining structure is identified with a principal 
fiber bundle with the structure group 5*0(3). The orbits of this group are the fibers 
and the base space is what we have called the reduced configuration space (shape space 
in the language of |^5| ) . Two kind of frames play an essential role in this approach: 
the space frame, which is an externally prescribed inertial frame, and the body frame, 
which can be attached to a flexible body in many ways. This arbitrariness in the choice 
of the body frame is a gauge freedom, any choice of gauge corresponding to a section 
of the fiber bundle. A gauge potential A a is introduced. The decomposition of the 
motion in rotational and vibrational part is gauge dependent, however the curvature 
form, 

ab ~ dq a dq b 3k A a A b , (M) 

called also Coriolis tensor is gauge covariant. Here {q a } are generic coordinates in 
shape space and the indices a, b symbolize the tensorial character of the Coriolis tensor 
with respect to changes in the shape space coordinates, while indices i, j, k denote the 
components in some body frame. Thus the Coriolis tensor is a geometric object with 
3 x (37V - 6) x (3iV - 6) components. 

A metric G a b in shape space is also introduced by means of the horizontal 
contribution to the kinetic energy. By a Kaluza-Klein type decomposition the Riemann 
tensor and curvature scalar associated to this metric are found, respectively: 

2 V " '■' J ~' , ' d ±J "-' '- , '! 1 "'" 7 ' 

3 

— 1 

4 

In the second formula the indices a, b were raised by the inverse of the (non-flat) metric 

G a b- 

Specializing to vanishing angular momentum leads to horizontal motions in the 
principal fiber bundle. It is thus not surprising, that we have exploited a somewhat 
simpler geometric structure in our analysis. 

Our approach differs in many aspects, as we were discussing a purely relational 
theory. We cannot start from the concept of space frame. It would imply to monitor 
the evolution of the system with respect to some preexisting inertial system, which is 
not the case. By contrast, we have chosen a definite body frame for our computations 
in Sec. 4. This is the principal axis frame. 

We have imposed constraints on the dynamics by means of Lagrange-multipliers, at 
the end of the reduction process arriving to a Jacobi- type action (|29|). The Lagrangian 
in |p5| , without the kinetic energy of the center of mass and specialized to vanishing 
angular momentum (when the vertical contribution to the kinetic energy about the 
center of mass vanishes) is 

L(q,q) = ^G ab q a q b -V(q) . (56) 



Rabcd = - [B i ab Ii j B 3 cd - B a[c I M B^ d]b j (54) 
R = lB i °»I ij Bi b . (55) 
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By writing the Jacobi principle version of this action with the energy set to zero we 
can identify the metric G a b with our metric GAaBb, expressed in a coordinate system 
adapted to the six generators. 



By comparing our formulae (49) and (BOJ) for the Riemann tensor and curvature 



scalar with Eqs. (54) and (55) respectively it is tempting to identify our vorticity 
tensor cj^sfe with the Coriolis tensor. For this purpose first we write Eq. (55) in the 
principal axis frame: 

/ T TDVab I T TDV 



Second, we remark that unlike the Coriolis tensor, the vorticity tensors are defined as 
three geometric object with 3N x 3iV components. However by their definition (42) as 
projections we see that the vorticity tensors live in the (3iV — 6)— dimensional reduced 
space (shape space). When coordinates adapted to this space are chosen, they will 
have equally (3 AT — 6) x (37V — 6) nonvanishing components. Third, indices Aa, Bb 
are raised with the (inverse) flat metric g AaBb ; in contrast with the contravariant 
form of the Coriolis tensor in (^7|), where indices are raised with the nonflat metric 
G ab . But there is no difference in raising the indices of ^>AaBb with the nonflat metric 
QAaBb e it ner _ This is because the two metrics differ only in terms homogeneous in 
the rigging vectors. When contracted with the vorticity tensor (a projected object) 
these difference terms vanish by virtue of Eqs. (f4l|). 

Thus we identify the set of the three vorticity tensors oj^ with \JTZB V /2 (no 
summation). The advantage in employing the vorticity tensors, which are defined on 
the (unreduced) configuration space is that their indices can be raised and lowered 
with the flat metric g. 



6. Concluding Remarks 

The study of the Barbour-Bertotti theory describing the Newtonian N— particle 
system constrained by the vanishing of all constants of the motion revealed that a 
complete reduction at the Lagrangian level is possible whenever the particles are not 
along a line. 

The reduced space emerged as the space of orbits of the Leibniz-group. Free 
motions are geodesies with respect to the Rlemannian metric GAaBb while motions 
characterized by V ^ const are geodesies of the conformally scaled metric 
— 2V{x)GAaBb- The vorticity tensors of the generators of rotations were introduced 
and their expression in terms of the pricipal moments of inertia was derived. 

The Riemann tensor and the curvature scalar associated with the metric GAaBb 
were given in terms of the vorticity tensors. The additional expression of the scalar 
curvature in terms of the principal moments of inertia allowed for a study of the 
collinear configurations. Curvature singularities arise in these configurations, unless 
the number of particles is three. 

The Riemannian curvature tensor and the curvature scalar associated to the 
conformally scaled metric GAaBb = —2V(x)GAaBb, which applies to the case V ^ 
const can be found by applying the relations among these objects for conformally 



related metrics [30 . For illustration we give here the curvature scalar R associated to 



14 



the metric G: 

V A aV Bb log(-V) 



R = — \ R - (3N - 7)G AaBb 



-IV 



3N 

— - 2 ) V Aa \og{~V)V Bb log(-F) 



(58) 



4 

For any definite choice of the potential V, such relations can be further exploited. 
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